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^ [ Abstract 
vC 

We prove that the maximum number N c of non-relativistic electrons that a nucleus 
of charge Z can bind is less than 1.22Z + 3Z 1 / 3 . This improves Lieb's upper bound 
N c < 2Z + 1 [Phys. Rev. A 29, 3018-3028 (1984)] when Z > 6. Our method also 
applies to non-relativistic atoms in magnetic field and to pseudo-relativistic atoms. We 
^ ■ show that in these cases, under appropriate conditions, limsup^^^ N c /Z < 1.22. 
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■ 1 Introduction 

t> 

^ \ Let us consider an atom with a classical nucleus of charge Z and iV non-relativistic quantum 

electrons. The nucleus is fixed at the origin and the iV-electron system is described by the 
Q\ ■ Hamiltonian 

i=l v 1 41 7 l<i<j<N 1 * 3 1 



O 



acting on the antisymmetric space f\ (L 2 (R 3 ) <8> C 2 ). The nuclear charge Z is allowed to be 

5-H ' 



1=1 

any positive number, although it is an integer in the physical case. 

The ground state energy of N electrons is the bottom of the spectrum of Hn,z, 

E(N, Z) = inf spec H N<Z = inf (i[>, H N)Z ip). 

I M 1.1,2=1 

We say that N electrons can be bound if E(N, Z) < E(N — 1,Z), namely one cannot remove 
any electron without paying some positive energy. Due to the HVZ theorem (see e.g. [27J, 
Theorem 11.2) which states that 

ess spec H NZ = [E(N — 1), oo), 

one always has E(N) < E(N — 1). Moreover the binding inequality E(N, Z) < E(N — 1, Z) 
means that E(N, Z) is an isolated eigenvalue of H N z . Zhislin (1960) [28] show that binding 
occurs provided that N < Z + 1. 
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Of our interest is the maximum number N c = N C (Z) of electrons that can be bound. It 
is a long standing open problem, sometimes referred to as the ionization conjecture (see e.g. 
PH [H [251 12]), that iV c < Z + 1, or maybe N c < Z + 2. Note that N c > Z due to Zhislin's 
result. We now briefly present the status of the conjecture, and we refer to [12] (Chap. 12) 
for a pedagogical introduction to this problem. 

It was first proved by Ruskai (1982) [IB] and Sigal (1982, 1984) [2D [22] that N c is not 
too large. In fact, Ruskai [IB] showed that N c = 0(Z 6 / 5 ) as Z — y oo and Sigal [22] showed 
that N c < 18Z and limsup^oo N c /Z < 2. Then Lieb (1984) [10] gave a very simple and 
elegant proof that N c < 2Z + 1 for all Z > 0. Lieb's upper bound settles the conjecture for 
hydrogen but it is around twice of the conjectured bound for large Z. 

For large atoms, the asymptotic neutrality lim^-^oo N c /Z = 1 was first proved by Lieb, 
Sigal, Simon and Thirring (1988) [1.1]. Later, it was improved to N c < Z + 0(Z 5 ^ 7 ) by 
Seco, Sigal and Solovej (1990) [19] and by Fefferman and Seco (1990) [7]. The bound N c < 
Z + const, for some Z-independent constant, is still unknown, although it holds true for 
some important approximation models such as Thomas- Fermi and related theories [Bl [3] and 
Hartree-Fock theory [24] [25]. 

In spite of the asymptotic neutrality, Lieb's upper bound N c < 2Z + 1 [TO] is still the 
best one for realistic atoms (corresponding to the range 1 < Z < 118 in the current periodic 
table). The purpose in this work is to find an improved bound for all Z > 0. As in [10], we 
do not need the binding inequality; more precisely, that E(N, Z) is an eigenvalue of H^ z is 
sufficient for our analysis. One of our main result is the following. 

Theorem 1 (Bound on maximum ionization of non-relativistic atoms). Let Z > (not 
necessarily an integer). If E(N,Z) is an eigenvalue of H^ t z, then either N = 1 or 

N < X.22Z + 3Z 1/3 . 

The factor 1.22 can be replaced by ft^ 1 with (3 being defined by (T^j below. 

Remark 1. The bound 1.22Z + 3Z 1 / 3 is less than Lieb's bound 2Z + 1 when Z > 6. 

Remark 2. While Lieb 's result holds true for both fermions and bosons, our result only holds 
for fermions (in fact, our method works also for the bosonic case but it yields an estimate 
worse than Lieb's). Note that the ionization conjecture only concerns fermions since for 
bosonic atoms it was shown that lim^^oo N c / Z = t c w 1.21 by Benguria and Lieb [2] and 
Solovej I23tf (the numerical value 1.21 is taken from JT$). In our proof below, we use Pauli's 
exclusion principle in Lemma® More precisely, we use the fact that in a fermionic atom the 
average distance from the electrons to the nucleus of charge Z is (at least) of order Z^ 1 ^ 3 . 
In contrast, the corresponding distance in the bosonic atoms is of order Z~ x . 

Remark 3. Although Lieb's method Ffflf can be generalized to molecules, we have not yet 
been able to adapt our method to this case. 

Our method also applies to other models such as non-relativistic atoms in magnetic 
fields and relativistic atoms, and we shall discuss these extensions later. In the rest of the 
introduction let us outline the proof of Theorem [I] As a first step we get the following 
bound. 
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Lemma 1. If E(N, Z) is an eigenvalue of H^,z then we have 



a N (N - 1) < Z(l + 0.68 N- 2/3 ), 

where 

a 7v := mi 



x 1 ,...,x N e<">- N 



(jv-i)e W 



i=l 



This result is proved by modifying Lieb's proof: in [10] Lieb multiplied the eigenvalue 
equation (Hn^z — E(N, Z))^n,z — by \xn\^n,z- We instead multiply by x 2 N ^^^z and 
employ the Lieb-Thirring inequality to control error terms. 

Roughly speaking, the number a^ 1 yields an upper bound on N/Z. This bound improves 
previous results since is bigger than 1/2 (one can see that a 2 = 1/2 and qjjv > a/5/4 = 
0.559... when iV > 3). Although we do not know the exact value of a>N, it is possible to 
derive some effective estimates. We may think of ajv as the lowest energy of N classical 
particles acting on IR 3 via the potential V(x, y) = t~T; under some normalizing condition. 
It is natural to believe that if N becomes large, then converges to the statistical limit 

S d p( x ) My) 




2\x-y 



(3 := inf < „ : p a probability measure on R > . (2) 

J \x\ ap(x) 

Results of this form in bounded domain have already appeared in [T7j. Indeed, we can show 
that ct/v actually converges to /3 and provide an explicit estimate on the convergence rate. 
Theorem [1] essentially follows by inserting the lower bound on in Proposition [1] below 
into the inequality in Lemma [TJ 

Proposition 1. The sequence {«at}/v=2 ^ s increasing and for any N > 2 we have 

(3>a N > J ^[(3-3((3/6) 1 / 3 N^% 
with (3 being defined by (TJj. Moreover, f3 e [0.8218,0.8705). 

Remark 4. We do not know the exact numerical value of (3, but our bound that (3 G 
[0.82188, 0.8705) is already rather precise. There is of course still room for improvement. 

The article is organized as follows. We shall prove Theorem [1] in Section 2. Then we 
discuss some possible extensions of our method in Section 3. Proposition [1] is of independent 
interest and we defer its proof to Section 4. 



Acknowledgments. I am grateful to Mathieu Lewin and Jan Philip Solovej for very helpful 
discussions, and to Rupert L. Frank and Elliott H. Lieb for pointing out the lower bound 
( flQj) which improves the constants in Lemma |2] and Theorem [lj I thank the referee for his 
constructive comments. This work was done when I was a visiting student at Departement 
de Mathematiques, Universite de Cergy-Pontoise (France), and I wish to thank the people 
there for the warm hospitality. 
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2 Proof of Theorem [D the new bound 



2.1 Lieb's method 

In order to make our argument transparent we start by quickly recalling the proof of Lieb 
[TO] . Assume that E(N,Z) is an eigenvalue of Hn,z corresponding to some normalized 
eigenf unction ty N . Multiplying the Schrodinger equation 

{H Ni z-E(N,Z))* NtZ = (3) 

by \xn\^n,z and then integrating, one gets 

= {\x N \* N ,z,(H NtZ -E(N,Z))* NtZ ) 

= (\xn\^n,z, {Hn-i, Z - E(N, Z))^ff iZ ) + 2 {\xn\*n,z, -A n V n , z ) 



N.Zi 



N-l I | 
\Xn\ 



1=1 



\Xi — Xn\ 



^n,z ) ■ (4) 



The first term in the right hand side of (|1J) is non-negative since i/jv-i.z > E(N — 1, Z) > 
E(N, Z) (in the space of iV — 1 particles x\ v ..,Xff-x)- The second term is also non-negative 
due to the inequality 

Re (\x\f, -Af) > for all / G H^M 3 ). (5) 

Thus the third term in (j4]) must be non-positive. Using the symmetry of I^jv^I 2 (which 
holds true for both fermions and bosons) we can rewrite it as 



l<i<j<N 1 1 jl 



I 7 



It follows from the triangle inequality that 



1 y- N±M > 1_ (6) 

l<i<j<N 1 n 

Hence we obtain —ZN + Af ( J ^~ 1 ) < o, namely iV < 2Z + 1. The inequality is strict since the 
triangle inequality is strict almost everywhere in (M 3 )^. Note that the lower bound 1/2 in 
(Q is sharp (when \xi\ <C \xj\ if i < j). 

2.2 Proof of Lemma [I] 

Instead of multiplying the equation (EJ) by \xn\^n,z, we now multiply by x 2 n ^n,z and inte- 
grate. We obtain 

= (x 2 N V N;Z , (H N _ 1;Z - E NiZ )^ N)Z ) + - (x 2 N V N;Z , -A N V NiZ ) 



+ U NiZ l-z\x N \ + ± 

\ V l<i<i5 



Xf + Xj 



i<j<N 1 * 



N,Z 



> o (^iv.z, -A N V NtZ ) + (*ar,z, (-^ + «jv(N - 1)) |xjv|*jv,z) • 



Recall that is defined in ([T]). This implies that 

a N (N -l)<Z-~ (x 2 N ^ N>z , -A N V N>Z ) (*n,z, M^n^)' 1 . (7) 

As we will see, the main advantage of our method is that the number a at is bigger than 
1/2 when N >3. However, we do not have an inequality similar to (jSJ) with \x\ replaced by 
x 2 . In fact, for all / G ff 1 (IR 3 ), applying the identity 



Re { V f, -A/) = ( yV 2 f, ( _ A 
to (f(x) = x 2 we find that 



Re<x 2 /,-A/> = (f,(\x\(-A)\x\-l)f) > ~(f,f) (9) 

by Hardy's inequality, —3/4 being the sharp constant. 

Our observation is that we may still control the second term in the right hand side of (jTJ) 
since ($n,z, \xn\^ tv,z) _1 is small (in comparison with Z) . In fact, ($n,z, \xn\^n,z) can be 
understood as the average distance from N electrons to the nucleus, which is well-known to 
be (at least) of order Z~ l l 3 . We have the following explicit bound. 

Lemma 2. If^>N,z is a ground state of H^ z then 

(9 N>Z , \x N \^ N>z ) > 0.553 Z~ l N 2/3 . 
It follows from (Q and Lemma [2] that 

~ (a? N * N ,z, -M> N , Z ) (* N;Z , IxnI^^z)- 1 > -0.68 ZN^ 3 . 

Substituting the latter estimate into © we obtain the inequality in Lemma [TJ We now 
provide the 

Proof of Lemma\j^ The following proof essentially follows from [12] ( p. 132). Note that 

(®N,z, \x n \^n,z) = 77 / \x\py NZ (x)dx 
where the density p^ NZ of ^n,z is defined by 

P*n,z( x ) '■= N ••• / \^N,z(x,o-i;x 2 ,a 2 ; ...;x Nj a N )\ 2 dx 2 ...dx 



(71=1,2 a N = l,2 E 3 (JV _ 1) 



By solving the Bohr atom as in (9] (after eq. (40) p. 560) one has the lower bound on 
the ground state energy 



(10) 



where A = (3 1/3 /2)2 2 / 3 . Moreover, one has the Lieb-Thirring kinetic energy inequality 

i - r 

: = 2 E - A ^^v^) ^/ P^, z (x) 5/3 rfx, (11) 



where if = 2" 2 / 3 (3/10) (2/(5L)) 2/3 with L = (vrS 3 /^)" 1 = 0.01225... (this constant L is 
taken from [B]). Since £?(iV, Z) = —)&& NZ by the Virial Theorem, we get from ( ITU]) and ( ITTj) 
that 

J P* N , z ( x ) 5/3dx ^ K-'AZ 2 ^ 3 . (12) 
On the other hand, we have the following inequality introduced by Lieb ([9], p. 563) 

p/2 / \ / \ l+5p/6 

y?(x) 5/3 dx I f / |x|V(z)dx > C p f / ^(x)dx 



for any nonnegative measurable function (f(x), with the sharp constant C p being attained 
for <p{x) = (1 — |a;| p )+ 2 . In particular, applying this inequality to (p(x) = p^, N z (x) and p — 1, 
we get 

1/2 

p* JVZ (a;) 5/3 dx I / |z|p^(x)da; > diV 11/6 (13) 



;>:{ 



where Ci = 7r -i/3 2 -i35/3 5 5/6 7 i/3 11 -3/2 = .4271... Combining (Q2D and ([33]) we obtain the 
desired inequality. □ 

2.3 Proof of Theorem Q] 

Let us admit Proposition [ljfor the moment and derive Theorem [TJ Lemma [1] and Proposition 
UJ together yield a lower bound on Z in terms of N, 

N(/3 - 3((3/6)^N-^) 

1 + 0.68 iV-2/3 < Z - (14) 

It is just an elementary calculation to translate ( !T4|) into an upper bound on N in terms of 
Z. If mm{N,Z} < 3, then max{2, fi^Z + 3Z 1 ^} >2Z + 1 (since < 0.8705), and hence our 
bound follows from Lieb's bound N < 2Z+1. Ifmin{iV,Z} > 3, then pS^/Qy^N- 2 / 3 > 
and Lieb's bound implies that N/Z < 2 + Z^ 1 < 7/3. Thus the desired result follows from 
(JUJ) and the following technical lemma whose proof is provided in the Appendix. 



Lemma 3. For Z > 0, N > 0, N/Z < 7/3 and (3 > 0.8218 one has 

i w, f 1 + 0.68A^~ 2 / 3 I 

r ^ + 3Z / 3 >min |^___j. 
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3 Some possible extensions 
3.1 Atoms in magnetic fields 

In this section, we consider the ionization problem with the presence of a magnetic field. 
The system is now described by the Hamiltonian 

Z\ ^ 1 



i=l x 1 l|/ l<i<j<N 



| X i X j | 



N 

acting on the fermionic space /\(L 2 (M 3 ) eg) C 2 ). The kinetic operator is the Pauli operator 

T A =\a- (-iV + -4(a;))| 2 = (-iV + -4(x)) 2 + a ■ B, 

where A is the magnetic potential, B = curl (.4.) is the magnetic field and a = (a 1 , a 2 , a 3 
are the Pauli matrices 

i ( 1 \ 2 f -i \ 3 ( 1 

For simplicity we shall always assume that A G Lf oc (IR 3 , R 3 ), V-i £ L 2 oc (R 3 ) and 
|£>| G L 3//2 (R 3 ) + L°°(R 3 ). Under these assumptions, it is well known that (— iV + ^4(s)) 2 is 
essentially self-adjoint on L 2 (R 3 ) with the core C^°(IR 3 ) [TB], and \B\ J r Z /\x\ is infinitesimally 
bounded with respect to (— iV + A(x)) 2 (see e.g. [2D])- In particular, the ground state energy 

E(N, Z, B) = inf spec H N>Z>A 

is finite. We shall also assume that iV i-> J5(iV, Z, B) is non-increasing (for example, this is 
the case if B = (0, 0, B) is a constant magnetic field [33]). Note that the ground state energy 
depends on A only through B by gauge invariance (see e.g. [12J p. 21). 

Of our interest is the maximum number N c such that E(N C , Z, B) is an eigenvalue of 
Hn,z,a- Seiringer [20] showed in 2001 that 

Nc<2Z+1 + ^,Z,B)-E(N kZ,B) 

2 N c Z(k-l) V ; 

for all k > 1. In the homogeneous case, B = (0, 0, B), his bound yields 

N c < 2Z + 1 + CiZ 1/3 + C 2 Zmin{(B/Z 3 ) 2/5 , 1 + | ln(£/Z 3 )| 2 } . (16) 

In particular, in the semiclassical regime lim^_ s>0O (5/Z 3 ) = 0, Seiringer's bound implies that 

N c 

limsup — < 2. 

Z->oo Z 

In contrast, it was shown by Lieb, Solovej and Yngvason (1994) [13] that if lim^_ s . 00 (_B/Z 3 ) = 
oo, then 

liminf > 2. 

We shall improve these upper bounds using the method in the previous section. Our 
result in this section is as follows. 
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Theorem 2 (Bounds on maximum ionization of atoms in magnetic fields). Assume that 
A E Lf 0C (M 3 ,M 3 ), V • A e L 2 0C (R 3 ) and \B\ = |curl(.A)| e L 3/2 (M 3 ) + L°°(M 3 ). For any 
Z > 0, denote by N c = N C (Z) the maximum number such that E(N C , Z, B) is an eigenvalue 
of Hn,z,a- Then we have, for every k > 1, 



(17) 



If B = (0,0, B) is a constant magnetic field, then 

N c < (1.22 Z + 3Z 1/3 ) (l + 11.8 Z- 2/3 + min jo.42 (B/Z 3 ) 2/5 , C(l + | ln(5/Z 3 )| 2 )} 

/or some universal constant C (independent of Z and B). In particular, if \\mz->oo{.B / Z 3 ) = 
0, then 

limsup-^ < 1.22. 

Z^>oo Z 

The number 1.22 in all bounds can be replaced by with (3 being defined by |Ij). 

Proof. Assume that ^n,z,a is a ground state of B nz ,a- Following the proof of Lemma [TJ we 
have 

a N (N -1)<Z- (x 2 N t! N>Z)A , T a V NiZ ,a) (^n,z,a, \xn\^n,z,a)~ X > ( 18 ) 

which is the analogue of (JTj). 

We may assume that N > @~ l Z + 3Z~ 2/3 (otherwise we are done). In this case the left 
hand side of ( ITS]) can be bound by 

*MN-l)> g _ 1+ N 3z _ 2/r (19) 

This estimate follows from the lower bound on oln in Proposition[T]and the following technical 
lemma whose proof is provided in the Appendix. 

Lemma 4. For Z > 0, N <E N, N > p~ x Z + 3Z~ 2 / 3 and f3 > 0.8218, one has 

{P - 3(/3/6) 1/3 iV- 2/3 )(/3- 1 + 3Z- 2 ' 3 ) > 1. 

The second term in the right hand side of ([TBI can be bound in the same way as in [20J. 
More precisely, using (jSJ) with —A replaced by T A > 0, one has 

(x 2 n ^n,z,a^ t a^n,z,a) = (*n,z,a, (\xn\Ta\xn\ - 1) *n,z,a) > -1- (20) 
On the other hand, for every k > 1, 

(^N,Z,A, IxnI^N^a)' 1 < (^N,Z,aA x n\~ 1 ^N,Z,a) 

_ i^N,Z,A' Hn,Z,A^ 'n,Z,a) — (^N,Z,A' Hn,UZ,A^ N,Z,a) 



NZ(k- 1) 
E(N. Z. B) - E(N. kZ. B) 



since ^n,z,a is a ground state of H nz ,a- Then (JTTJ) follows by substituting ffT^l) . fl2"Ul) and 
d2U) into (OS]). 

Now assume that B = (0, 0, B) is a constant magnetic field. It follows from [T4] (Theorems 
2.4, 2.5) that if iV > Z/2, then the ground state energy E(N,Z,B) := E(N,Z,B) can be 
bounded from below by 

E(N,Z,B) > -NZ 2 (l8.7Z" 2/3 + min jo.95 (B/Z 3 ) 2/5 , C (l + | \n(B/Z 3 )\ 2 ) }) (22) 

for some universal constant C (independent of N, Z and I?). (It is obtained when applying 
(2.27), (2.26), (2.29) in p3] to the cases: B < Z 4 / 3 , B > Z A ''\ B > Z 3 , respectively.) 

We can choose k = 2 in (IT7j) . Then the desired bound follows by using the upper bound 
E(N, Z, B) < and the lower bound on E(N, 2Z, B) derived from (j22D- □ 

AT 

Remark 5. We may also consider the Hamiltonian H^z.a on the bosonic space <g) (L 2 (R 3 )® 

sym 

C q ), where q is a spin number. In this case the inequality fn\ ) still holds true. Moreover, if 
B = (0,0,5) is a constant magnetic field, then using the estimate (p. 1948) 

E(N,Z,B) = NZ 2 E(1,1,B/Z 2 ) > ~NZ 2 min {l + AB/Z 2 , C\ \xi(B/Z 2 )\ 2 } 



we get from (fI7P that 

N c < {P^Z + 3Z 1 /' 3 ) (1 + min {1 + 4B/Z 2 , C 2 \ \n(B/Z 2 )\ 2 }) . 
In particular, if \im z ^oo(B / Z 2 ) = 0, then our bound yields 

limsup— < 2/T 1 < 2.44. 

Z^oo Z 

It slightly improves the bosonic bound in f2(% . which gives limsup z _ >00 (A'" c /Z) < 2.5. 



3.2 Pseudo-relativistic atoms 

In this section we consider the pseudo-relativistic Hamiltonian 

< z = E(«- 1 (^T^-a-')- 1 ||)+ £ ,-4-, 

N 

acting on the fermionic space /\(L 2 (IR 3 ) <g> C 2 ). Here a > is the fine-structure constant. It 
is well known that the ground state energy E iel (N, Z) := inf spec Hffi z is finite if and only if 
Za < 2/n (see e.g. [12]). The physical value is a = e 2 /(he) ~ 1/137 and hence Z < 87.22. 
However, we allow a to be any positive number. 

As in the previous dicussions, we are also interested in the maximum number N c such 
that the ground state energy E Tel (N c , Z) is an eigenvalue of H T ^ l cZ . Note that Lieb's bound 
iV c < 2Z + 1 still holds in this case. In fact, due to a technical gap the original proof of Lieb 
in [TD] works properly only when Za < 1/2. However, it is posible to fill this gap to obtain 
the bound up to Za < 2 /tt [5J. On the other hand, to our knowledge, no result about the 
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asymptotic behavior of N c /Z is available for the pseudo-relativistic model, although within 
pseudo-relativistic Hartree-Fock theory it was recently shown by Dall'Acqua and Solovej 
(2010) [U that N c < Z + const. 

Our result in this section is the following. 

Theorem 3 (Bound on maximum ionization of pseudo-relativistic atoms). Let Z > such 
that Za < k < 2/tt. If E rel (N, Z) is an eigenvalue of H r ^ x z , then either N = 1 or 

N < 1.22Z + C K Z 1/3 

for some constant C K depending only on k. The number 1.22 can be replaced by with f3 
being defined by (TJ|). 

Proof. Assume that ^ r ^ z ls a ground state of H r ^ l z . As an analogue of (|7|) we get 



a N (N - 1) < Z - (x 2 N W$ z , crV-Av + cr 



The left hand side of (1231) can be bounded using ([191) . Turning to the right hand side of 
231) . we first show that for any function / : M 3 — > C smooth enough 



Re ( x ,a 



1 (v^A 



a: 



-i 



a 



-i 



./ 



L 2 (R 3 ,dx) 



> -o </»/>■ 



(24) 



It suffices to show ( 1241) for a = 1 (the general case follows by scaling). Using the Fourier 
transform f{p):= J e~ t2wp ' x f (x)dx and applying §8$) to 



we find that 



<p(jp) := Vi^P) 2 + 1-1 



(27T P f 



Re(x 2 f, 



-A„, + 1-1 



./ 



^{2ixp) 2 + 1 + 1 



L 2 (R 3 ,dx) 



(27r)- 2 Re(-A p /,^/ 
(2vr)- 2 ( ^ 2 f, (-A, 



L 2 (R 3 ,dp) 



2lp 



V 1/2 f 



Ap 2 ((2np) 2 + 1 

Then it follows from Hardy's inequality — A p > l/(4p 2 ) that 



Re(x 2 /, V-A* + 1-1 / 

J J I L 2 (R 3 ,dx) 

2 v /(2irp) 2 + 1 + 1 



> - / 



4((27rp) 2 + l)( v /(2 7 rp) 2 + 1 + 1 



■./' 



L 2 (R 3 ,dp) 



> 4 (it) =-l(fj) ■ 
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The term (fy,\x]\r\ty) 1 can be estimated similarly to ( 12T|) . namely 



for every fc > 1 such that £;Zct < 2/vr. It is well known that > E(N, Z) > -C K Z 7 /' 3 
provided that Za < k. In fact, it was shown by S0rensen [26] that, in the limit Z — > oo 
(and Za = k fixed), the leading order of the ground-state energy E(N, Z) is given by the 
Thomas- Fermi theory which is of order Z 7 ! 3 . Thus we can conclude that 

(^M^r 1 <C K Z- 2 '*. (25) 
The desired result follows from (121]), (j25) . p5]) and (115) . □ 



4 Proof of Proposition [D Analysis of ajy 

This section is devoted to the proof of Proposition [TJ For the reader's convenience, we split 
the proof into several steps. Recall that a at and (3 are defined in (JTJ and (j2J), respectively. 

Step 1. The sequence is increasing in N and it converges to (3 as N — )■ oo. 

Proof. The fact that a at is increasing is shown as follows: for every x\,...,Xn G M 3 we have 

+ x\ ^ f 1 ^ x\ - 



l X i X j\ V 2) |Xi Xj, 



l<j<j<Af 1 Jl k=l \ y ' i<j;i^k,j^k 




N 



n v-i(-/V — 1) |a- 

i=l 



where we have used the definition of a^-i- This implies that a^-i < oin- 

We shall show that converges to (3. We start with the upper bound < (3. Let p 
be an arbitrary probability measure on M 3 . Then 

ff x 2 + y 2 , . , . . , /" 1 z? + 

JJ 2W^V\ Mx)My) = J i^v^^R^^-.-dp^) 

R 3 xR 3 R 3iV 1S*<?SJV 

> y ^ ( W J dp(xi)...dp(xjv) = «tv y |x|dp(x). 

R 3iV \ i=1 / R 3 

Thus a N < (3 for all N >2. 

Let us prove a lower bound. For the reader's convenience, we give now a simple bound 
which is enough to get that converges to (3. We will provide a better lower bound in the 
next step. 

Let {xi}^ be N arbitrary distinct points in R 3 and let r > 0. For our purpose we 
may assume that 5^»=i \ x i\ = N. For every i, let dpi be the uniform measure on the sphere 
\x — X{\ = Ti with the radius Tj := r\xi\ such that J dpi = 1. Define the probability measure 

d P(X> := ^Eli d /^)- 
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Since J \x\dp(x) > 1 (due to the convexity J \x\dfjLi(x) > \xi\), we have 

x 2 dp(x) dp(y) 



\x - y\ 



On the other hand, 

x 2 dp(x) dp(y) 

\x - y\ 



hi 



x 2 dpi(x) dpj(y) 



< N- 2 (l + r) 2 J2 



x-y\ 

x 2 dpi(x) dp^y) 
\x - y\ 



< N- 2 (l+r) 2 



E 



rn 2 



CCj Ob q 



N 
+ — 



The first inequality follows from |a:| < (1 + r)\xi\ for every x on the sphere \x — Xi\ = r», and 
the second inequality is due to Newton's theorem (see, e.g. [12], p. 91). Thus 



This implies that 



E 



N 



x- 



OlJJ OC t 



> (1 +r)- 2 N 2 /3 -r^N. 



a N 



^^^[(l+r)- 2 /?-^)- 1 ] 



for all r > 0. 



(26) 



We can choose, for example, r = N 1 ^ to conclude that — > (3 as N — > oo. This ends the 
proof of Step 1. □ 

We now improve the lower bound (1261) . 

Step 2. We have the lower bound 

a N >JL l \p-3{p/6)^N- 2 ^] 
Proof. In fact, we shall prove that 



a N > 



N 



N-l 



1 + r 2 /3 1 



1 + r 2 







rN 



> 



N 



N-l 



1 



2r 2 



(27) 



for all r G (0, 1]. The desired result follows by choosing r = (4/3A^/3) -1 / 3 which maximizes 
the right hand side of ( 1271) . 

The bound ( |2"7|) is shown by following the same method as for (T26]) . but with more careful 
computations. We shall prove that (with the notation of the proof of StepHJ), for r e (0, 1], 



f r 2 
/ \x\dp(x) = 1 + — 



(2? 
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and 



x 2 dp(x) dp(y) 
\x - y\ 



< N~ 2 {l + r 2 



E 



x 2 



N 
+ — 



(29) 



The identity (|28p follows from a direct computation using the formula 

1 



f(x)dfii(x) 



\S 2 \ 



f(xi + riu)du 



s 2 

2tt 7T 



\S 2 \ 




f(xi + rj(cos 9, sin 9 cos (p, sin 9 sin p)) sm(9)d9 dp 



o o 



for any integrable function /. Here the second identity comes from the spherical coordinates 
u = (cos 9, sin 9 cos p, sin 9 sin p), where 9 G [0, 7r) and 9? G [0,27r]. (Note that if r > 1, then 
the left hand side of ([28]) becomes Nr(l + l/(3r 2 )).) 
Now we prove (12"9~1) . Using Newton's theorem we have 



x 2 dp(x) dp(y) 
\x - y\ 



x 2 dpi(x) d/Xj (?/) 



|x - y\ 



< ^- 2 E 

< A^~ 2 (l + r 2 



x 2 dpi(x) 



E 



Er^-r^ 

A¥=3 



2 N 



(1 + r 2 )x 2 dpi{x) 



Kl 



+ E ^ 



where 14: = and 



2xj.(x — Xi) d/ii(x) 



\S 2 



2Xj.TjU! 



\Xi - Xj + TiOol 



dui 



S2 



2 riXi(xi - Xj)mm{\xi - Xj\,rj} . f . ^ 
3 |xj — Xj|(max{|xj — rj}) 2 

Here we have used the formula 



\S 2 \ J \a + slo 



■duj 



s2 



1 a min{|a|,s} ~ 

'TTmT r xV fleR ' S>0 ' 

3 |a| (max{|a|, s\) z 



(30) 



Thus (|29|) will be validated if we can show that V^- + V,i > 0. We distinguish three cases. 
Case 1: \xi — Xj\ > max{r,, r 3 }. We have 



2 r 2 x i (x l -x j ) _ 2 r 2 x j (x j -x i ) 

3nr> . rp . 1 3 Q 'T* ■ T* ■ ^ 

r 2 (x 2 - x 2 ) 2 + (x 2 + x 2 )(xi - Xj) 2 



|3 



< 0. 
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Case 2: \xi — Xj\ < min{r^,rj}. In this case 

3 3 T j 3 T \ | X>i 1 . 1 30 j | J 

Case 3: < [rrr« — ar^j < r 3 - (the case Tj < \xi — Xj\ < is similar). We have 

2/y»2 ~, _ / , _ />-» _ \ O , >-» _ / , > ■ _ ry . ^ 



Vij + Vji 



I rp . <-y> .13 *3 /vi . 



It is obvious that < since \xj\ > \xi\. If V,- t j < then we are done; if Vy > 0, then using 

^ ^ Xj <^ Xj sj J 



r i < F« - Xj\ we get 



3 n 

It turns out that Vij + Vji > as in Case 2. □ 

We now turn to direct bounds on 0. 
Step 3. We have the bound 0.8218 < (3 < 0.8705. 

Proof. The lower bound follows from the following estimate whose proof will be provided 
later. 

Lemma 5. For any positive measure p on M 3 we have 

x 2 -\- y 2 

r dp{x) dp{y) 

\x - y\ 



> max | J j ^max{|x|, \y\} + ^ I » \ dp ( x ) djP ( y ) 



(min{ x 


■ 1 


y|}) 2 


max{ x 


1 j 


|y|} 



a? - 2/1 + o r I I i-| dp(x) dp(y) 



Remark 6. If p is radially symmetric, then three terms in LemmalE are equal. 

It follows from Lemma |5] that for any positive measure p on M 3 and for any A G [0, 1] we 
have 

x 2 + y 2 ff 

-dp(x)dp{y) > W x {x,y)dp(x)dp(y) 

\x-y\ JJ 

M 3 xR 3 R 3 xR 3 

where 

W x {x,y) := A ( max{|x|, \y\} + 



(min 


■{M 


•) \ 


y\}f 




x — 


y 





(min{ \x\ 


i \ 


y\}? 


max{ \x 


1 1 


\y\} 



(31) 
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It turns out that 

a ^ ■ f Wx(x,y) 
P > sup ml - — ; — -. 

As[o,i] ^ij/GK 3 \x\ + \y\ 

Thus the lower bound on /3 follows from the following lemma whose proof is provided in the 
Appendix. 

Lemma 6. With W\ being defined in l[3l\) one has 

sup inf ^4 > 0.8218. 

Ae[o,ip,yeR 3 \x\ + \y\ 

(A numerical computation shows that left-hand side of the inequality in Lemma [6] is equal 
to 0.8218066...) 

The upper bound on (3 is attained by choosing some explicit trial measure p. By restrict- 
ing p to radially symmetric measures we have 



oo oo 



r 2 dm(r)dm(s) 
max{r,s} 



if/ 

< A-ad := inf ^ — : m a probability measure on [0, oo) 

J rdm(i 



Choosing m{r) = fr 3 ' 2 l[i t g](r)dr, with dr being the Lebesgue measure, we get 

£ad<^-^ln(3) = 0.8704... 

(A numerical computation shows that /3 rat } is approximately 0.8702.) □ 

For completeness, we prove Lemma [51 

Proof of Lemma We start by proving 

x 2 + y 2 
\x - y\ 



(min 




i \ 


y\}? 




x — 


y 





dp(x) dp(y) > / / max{|z|, \y\} H , —. dp(x) dp(y). (32) 



We first show that (1321) follows from the following inequality: for any e > 0, if N is large 
enough, then 



Ei + 1 , ^ / r, , , ,-. minllxjl, \xA}\ , 

^—^->(l-e) J2 max{|x t |,|x J |}+ J.'.L.i ) (33) 

l<i<j<N 1 Jl l<i<j<N v 

for every C M 3 . In fact, we may assume that p(R 3 ) = 1. For every e > 0, taking N 
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large enough and using f[3"3"j) one has 



dp(x)dp(y) 



x 2 + y 2 



\x - y\ 

2 { E 



<X- ^ I j 



) 



N(N-1) \ " la* -a;,- 1 / 



dp(xi)... dp(xiv) 



> 



2(1 — e) v-^ / min{|xj|, |x 3 -|}\ 

N ( N _ ^ 2^ ^max{|xi|, \Xj\} + — \ x ._ x .\ — J Wi)- M^n) 



^ ^max{|x|, + 



min{ x 




x — 


y| 



dp(ar) dp(y). 



Since the latter inequality holds for every s > 0, the inequality (1521) follows. 

Now we show ( 1331) . This inequality follows from a key result of [H] . It was shown in [TT] 
(Theorem 3.1) that, for any e > 0, if iV is large enough, then 



max < V - N (} ^ J> > 

i<i<iv \xi — Xj\ \xA 



(34) 



for any {xi\f =l C M 3 . Since 
max{|xi|, \xj\} + 



(min 




) 


X 


,l}) 2 




x^ 


Xj 





< min{|xi|, \xj\} 



(max{ \xi 


i 


^■l}) 2 




X{ 


X 


3 





we can deduce from f l34"l) that 

— | - (1 - e) I max{|xj|, |xj|} + 



max < > 

1+3 



x 2 + x 2 



I 3/ 2 *^ j | 



(min 






X 


3\}? 




Xj 


Xj 





> 0. (35) 



Now take 1 > e > 0. For iV large enough, employing fl35|) repeatedly, we can assume that 



E 

l<i<i 



I ,y»2 

| X j I 



(1 — e) I max{|xj|, |xj|} + 



(min 




) 


X 


;l}) 2 




Xi 


Xj 





> 



for every eN < j < N. It turns out that 

c 2 ( 

— — (1 — e) I max{|xj|, |xj|} + 



E 

l<i<j<N 



X 2 + X 2 



(min{ Xj 


5 


X 


,l}) 2 




Xi 


Xj 





^ E 

l<i<j<eN 



X 2 + X 2 

X j X n 



(1 — e) ( max{|xj|, |xj|} + 



(min 




5 


X 


3\}f 




Xi 


Xj 





> — max{|xj|, |xj|} > — eN |xj| > — 



X<i<j<eN 



l<i<eN 



l<i<eN<j<N 



a* | , | Xj | } -|- 



l<i<j<iV 



(min 


{\Xi 


f 


X 


,l}) 2 




X{ 


Xj 
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Thus f[3"3"|) . and hence f[3"2"|) . follows. 

Next, we show that 

R 3 xR 3 R 3 xR 3 

This is equivalent to 

X ■ V dp W dp W > i // ^MMl Mx) My) . (37) 



R 3 xR 3 



\x — y\ 3 J J max{\x\,\y\} 



In fact, if p is radially symmetric, then fl36|) becomes an equality due to fl30|) . In the general 
case, let us introduce the positive, radially symmetric measure 



»*) = / 



SO(3) 



dR being the normalized Haar measure on the rotation group SO(3). Because of the positive- 

definiteness of the operator with the kernel , , , we can employ the convexity to get 

\ x y\ 



dp{x) dp{y) > / / X V dp{x) dp(y) 



\x — y\ J J \x — y 

R 3 x 

1 

3 

R^x 
1 



(min{|x| 


• 1 


y\}) 2 


max{ \x 


1 ? 


\y\} 



3 

R 3 xR 3 



(min{ \x\ 


•> \ 


y\}) 2 


max{ \x 


1 • 


\y\} 



dp(x) dp(y) 
dp(x) dp(y). 



Thus (1371) (and hence (136]) ) holds for all positive measures p. □ 



Appendix: Technical lemmas 

In this appendix we provide the proofs of some technical lemmas. 

Proof of Lemma\^ Let us denote f3\ := 3(/3/6) 1//3 for short. If the desired inequality fails, 
then /T 1 < N/Z < 7/3 and 

< 1 + 0.68 N~ 2/3 - (/T 1 + 3Z- 2/3 )(/3 - ^N- 2 ^) 
+ /3iiV- 2/3 (iV/Z - /T 1 - 3Z- 2/3 ) 
= N~ 2/3 [0.68 - 3f3{N/Z) 2/3 + Px{N/Z)] . 

Thus the polynomial 

h(x) := 0.68 - 3/3x 2 + /3 lX 3 

satisfies that h(-oo) = -oo, h(0) = 0.68 > 0, h^- 1 ' 3 ) < 0, h{{N/Zf/ 3 ) > 0, /i((7/3) 1 / 3 ) < 
and h(+oo) = +oo (to verify that ^(Z?" 1 / 3 ) < and /i((7/3) 1/3 ) < we need to use > 
0.8218). This implies that h(x) has more than three distinct roots, which is a contradiction. 

□ 
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Proof of Lemma\4\ Denote /9i := 3( ( 5/6) 1 / 3 . Assume that the desired inequality fails, namely 

3/3z -2/3 < ^-2/3^-1 + 3Z -2/3)_ (gg) 

Replacing the term + 3Z~ 2 / 3 in the right hand side of ( 138)) we get 

HI)- 

since /3 > 0.8218. Thus N > max{4Z, + 3Z~ 2 / 3 } > 4. 
On the other hand, (1381) is equivalent to 

N -2/3 > Z 3 1 



/?! 3/3(iV/Z) 2 / 3 ' 

Using /3 > 0.8218 and N/Z > 4 we have N < 4.5. It contradicts the fact that N must be an 
integer. □ 

Proof of Lemma® For any x,y G M 3 , denote a = max{\x\,\y\}, b = min{\x\,\y\} and 
c = \x — y\. Using the inequality u 2 + v 2 > 2uv for u, v > 0, we find that 

W x ( X ,y) = A (a + g + (l-A)(c + g 

= (A - A')a + f A'a + (1 - A)^ + (a^- + (1 - A)c^ 

> (A - A')a + f 2 y^A'(l - A) + 2 VA(1 - A) ) b 

for every < A' < A. We may choose A' such that 



A - A' = 2;/-A'(l - A) + 2y/X(l - A). (39) 



If A > 0.8 the solution to (BSD is 

; A + 2 



A' = 

Thus, for every x, y G M 3 , 



2v^T^A)-^(l-A)j 



The desired lower bound comes from g(0.843) = 0.821804... (A numerical computation shows 
that g(X) has a unique maximum at A = 0.843476... and g max = 0.8218066...). □ 
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